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Abstract 

The Hamiltonian structures of the incompressible ideal fluid, including entropy 
advection, and magnetohydrodynamics are investigated by making use of Dirac's 
theory of constrained Hamiltonian systems. A Dirac bracket for these systems is 
constructed by assuming a primary constraint of constant density. The resulting 
bracket is seen to naturally project onto solenoidal velocity fields. 

1. Introduction 

From the early work of Lagrange [1] it became clear that ideal fluid systems 
possess the canonical Hamiltonian form when one adopts a fluid element descrip- 
tion, the so-called Lagrangian variable description. Because the Lagrangian 
description is particle-like in nature, it is amenable to action functional and 
Hamiltonian formulations. However, when Eulerian variables are incorporated 
the canonical Hamiltonian structure for all ideal kinetic and fluid theories is 
altered because the transformation from Lagrangian to Eulerian variables is 
not canonical. This results in a Hamiltonian theory in terms of noncanonical 
Poisson brackets (see, e.g., [2, 3, 4, 5, 6, 7] for review). 
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The present paper concerns the proper treatment of the incompressibiUty 
constraint of fluid mechanics in the context of the Eulerian Hamiltonian the- 
ory in terms of noncanonical Poisson brackets. We do this by applying Dirac's 
method for incorporating constraints in Hamiltonian theories, a central element 
of which is a Dirac bracket. In the past, researchers have used Dirac brackets for 
various reasons in fluid mechanics [8, 9, 10, 11, 6, 12], but the flrst works to use it 
to explicitly enforce the incompressibility constraint for Euler's equation in three 
dimensions appear to be Refs. [13, 14, 15]. Here we first extend the work of these 
authors by constructing the Dirac bracket for the ideal fluid with the inclusion of 
entropy advection, which allows for the inclusion of any advected quantity like 
salt concentration in the ocean. This generalization reveals that Dirac brackets 
of the kind considered in Refs. [13, 14, 15], as well as our generalization, can 
be written in a considerably simplifled and perspicuous form in terms of the 
projection operator that takes a general vector fleld to a solenoidal one. With 
this realization we then construct the Dirac bracket for incompressible magne- 
tohydrodynamics (MHD), thereby making clear its Hamiltonian structure. We 
present these results together by starting from the full compressible ideal MHD 
equations, 

V = -V Vv - p-iy (p'^) + P"'(V X B) X B , (1) 
P = -V • {pv) , (2) 
B = V X (v X B) , (3) 
s = -V • Vs , (4) 

where v(x, t) is the velocity field, p(x, t) is the mass density, B(x, is the 
magnetic field, and s(x, i) is the entropy per unit mass. All of these dynamical 
variables are functions of x e Z-/ c M'^ as well as time. We suppose boundary con- 
ditions are such that no surface terms appear in subsequent calculations which, 
e.g., would be the case on a periodic box or all space. The observables of the 
MHD system are functional of these fields, denoted generically by F[p, v, B, s]. 
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In terms of these variables, this system has the following Hamiltonian (energy): 

H[p, V, B, s] = j Sx + pU{p, s) + ^) , (5) 

where = |vp and = |Bp. With the MHD noncanonical Poisson bracket 
of Refs. [16, 2] 

{F,G} = - j d^xiFpV ■G^ + F^-VGp-p-\'^ xv)-{F^xG^) 

+p-iVs • (F.G^ - F^Gs) + (p-'Fv • [VGb] - p-'Gv • [VFb]) • B 
+B . ([V (p-'F^)] • Gb - [V (p-'G^)] • Fb]) , (6) 

where F^ denotes the functional derivative of F with respect to v, i.e. F^ = 
5F/5v, and the same holds for Fg, Fb and Fp. Here the notation a • [M] ■ b = 
b- (a - [M]) is a scalar explicitly given by aiMijbj (with repeated indices summed) 
for any vectors a and b and any matrix (or dyad) [M]. The bracket (6) with 
Hamiltonian (5) gives the MHD equations (l)-(4) in the form F = {F,H}. 
(Assmning V • B = 0.) 

The paper is organized as follows: In Sec. 2 we review Dirac's formalism for 
constrained Hamiltonian systems. Then, in Sec. 3 this theory is used to ob- 
tain the noncanonical Poisson-Dirac bracket for the incompressible ideal MHD 
equations including entropy advection. Here we impose a primary constraint 
that is a constant and uniform density and the rest follows from Dirac's algo- 
rithm. In particular, it is seen that the corresponding secondary constraint is 
that the velocity field be solenoidal. We verify that Poisson-Dirac bracket in- 
deed produces the correct equations of motion. This is followed in Sec. 4 by a 
detailed comparison to previous attempts at incorporating incompressibility in 
Hamiltonian formulations of incompressible ideal fluids. Finally, in Sec. 5, we 
summarize and conclude. The paper also has several appendices that address 
various issues that arise in the text. 

2. Dirac brackets 

As stated above, Dirac's theory is used for the derivation of the Hamiltonian 
structure of Hamiltonian systems subjected to constraints. Dirac constructed 
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his theory in terms of canonical Poisson brackets and detailed expositions of his 
theory can be found in Refs. [17, 18, 19, 20, 5]. However, it is not difficult to 
show that his procedure also works for noncanonical Poisson brackets (cf., e.g., 
an Appendix of Ref. [10]). In this section, we recall a few basic facts about Dirac 
brackets in infinite dimensions in the context of noncanonical Poisson brackets. 

If we impose K local constraints $a(x) = for a = 1, . . . , if on a Hamil- 
tonian system with a Hamiltonian H and a Poisson bracket {•,•}, the Dirac 
bracket is obtained from the matrix C defined by the Poisson brackets between 
the constraints, 

C„^(x,x') = {$a(x),$Mx')}, 

where we note that Ca^(x, x') = — C/3a(x',x). If C has an inverse, then the 
Dirac bracket is defined as follows: 

{F,G}. = {F,G} ~ jd'xj dV {F,$„(x)}C-;(x,x'){$^(x'),G}, (7) 

where the coefficients G~^(x, x') satisfy 

J dVG-^i(x,x')G;3^(x',x") = J dVGa^(x,x')G^-^^(x',x") = 5c«^5(x-x"), 

which implies G~^(x,x') = -G^^(x',x). 

This procedure is effective only when the coefficients G~^(x, x') can be 
found. If C is not invertible, then one needs, in general, secondary constraints 
to determine the Dirac bracket. The secondary constraint is given by the 
consistency equation which states that ^i(x) = for the Hamiltonian H + 
J d^a;u(x)$i(x). This translates into 

J £x{<^>l{^),H}^ii^)^o, (8) 

for all functions n such that 

j rf^x^(x)G(x,x') = 0. 

Here the weak equality w stands for an equality on the manifold defined by 
$i(x) = 0. Equation (8) gives the expression which has to be satisfied by the 
secondary constraint. 
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3. Dirac bracket for ideal incompressible MHD 

To construct the Hamiltonian theory of ideal incompressible MHD, the first 
(primary) constraint is chosen to be a constant and uniform density po, i.e. 

$i(x) =p(x) -po- 

However, the Dirac procedure can be performed for the case of a nommiform 
background density (sec Appendix A). Given that Cii(x, x') = 0, at least one 
secondary constraint is needed. This secondary constraint, denoted $2(x), is 
given by {$i(x), H} = which leads us naturally to 

$2(x) = V • V. 

From the Poisson bracket (6), we compute the elements Ca^(x, x') as 

Cii(x,x') = 0, 

Ci2(x,x') = A5(x-x'), 

C2i(x,x') = -A<5(x-x'), 

C22(x,x') = V - (/9-i(Vxv) X V5(x-x')). 

From these expressions, we obtain the coefficients C~^(x, x') as 

C7iY(x,x') = A-iy • (p-'(V X v) X VA-i<5(x - x')) , 
Cr2^(x,x') = -A-i^(x-x'), 
C2Y(x,x') = A-i5(x-x'), 

C2-2^(X,X')=0, 

where A~^ acts on a function / as A~^/(x) = -(4?:)"^/ (i^x'/(x')/|x - x'|. 
Given the following expressions 

{$i(x),G} = -V-Gv, 

{$2(x), G} = ~AGp - V ■ {p-\\7 X v) X Gv) + V ■ (p-^VsG,) 
-V • (p-^VGb] • B) + V • (p-^V • [BGb]) , 
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we deduce various contributions to the Dirac bracket (7): 



J J rf3xdV{i?^,$i(x)}CiY(x,x'){$i(x'),G} = -y d3a;V-FvA-iV-(p-^(Vxv) X VA-V-C?v). 

yy d32.^3^'{F,$i(x)}Cf2^(x,x'){$2(x'),G} = j d^xV ■ (G^ + A'^V • (/9-^(V x v) x Gv) 
-A-^V ■ (p-^WsGs) + A-^V • (p-^VGb] • B) - A'^V • (p-^V • [BGb])) , 

J I rf3a;dV{i?^,$2(x)}G2-/(x,x'){$i(x'),G} = ' j '^'x (J> + A'^V • (p-^(V x v) x F^) 

-A-iy ■ (p-^VsF,) + A-iy • (p-^Vi^B] • B) - A-^V ■ (p-^V • [BFb])) V • Gv 

From the contributions associated with ^nd C21 , we notice that the part 
— / d^x{FpV ■ G^+F^ ■ VGp) of tlie Poisson bracket (6) vanishes. We also notice 
that the terms in the Dirac bracket only involve 

Gv := Gv - VA-^V • Gv =: P • Gv . (9) 

Two equivalent expressions for V acting on a vector a are P ■ a = a — VA~^ V • 
a = —V X (V x A~^a). The linear projection operator V acting on vectors is 
symmetrical, in the sense that 



j d^xa-V h = j d^xh-Va, 

for any vector fields a(x) and b(x). In addition, it satisfies the following prop- 
erties: 

V'^=V, V-V = 0, p.Vx=Vx, 

Vx7' = vx, v-r = o. 

As a consequence, we notice that the functional derivatives Gv are divergence- 
free, i.e. V-Gv = 0. In terms of Gv given by Eq. (9) the Dirac bracket is written 
in the following compact form: 

{F, G}, = J d^x (p-i(V X v) ■ (Fv X Gv) - p-'Vs • (F,Gv - F^G^) 

- {p-^F^ ■ [VGb] - p-^Gv • [Vi^B]) • B - B • ([V (p-^Fv)] • Gb - [V (p-^Gv)] • Fjip) 

Upon comparison with bracket (6), we see that this bracket is precisely that 
of Refs. [16, 2] with the functional derivatives F^ and Gv replaced by the 
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divergence- free functional derivatives and Gv according to Eq. (9). In this 
procedure, the terms of the bracket (6) in Fp or Gp disappear because V-Gv = 0. 
We also note that if we drop all terms but the first in Eq. (10), then with some 
manipulations one can show this bracket is equivalent to the one obtained in 
Ref. [13], albeit in a significantly simplified and perspicuous form, and that this 
term corresponds to the bracket of Ref. [21]. 

Because the Poisson bracket (10) is exactly the bracket of Ref. [16] with the 
replacement of the functional derivatives by projected functional derivatives, 
one wonders if one can always construct Dirac brackets by this procedure. In 
Appendix B it is shown that not all projections produce good brackets, only 
those that define Hamiltonian vector fields (see also Appendix C). 

Given that V • -FV = for all observables F, we obtain the following family 
of Casimir invariants of the Poisson bracket (10): 



where /(s) is any function of the entropy, i.e. it commutes with all the ob- 
servables, {C[s],G}* = for all G. This family originates from the family of 
Casimir invariants of the original Poisson bracket (6) given by / (fxpf{s), and 
the fact that the Dirac constraints are also Casimir invariants. This follows since 
Dirac brackets built on brackets with Casimir invariants retain those invariants 
(cf. Ref. [10]). As a consequence, the term J d^x pU{p,s) in the Hamiltonian 
is now a Casimir invariant, so that it can be dropped from the Hamiltonian 
because it will not give any contribution to the equations of motion (contrary 
to the compressible fluid or compressible MHD cases). Upon setting po = 1, the 
Hamiltonian becomes 



Therefore, the Hamiltonian theory for ideal Eulerian incompressible MHD is 
given by the bracket (10) with the Hamiltonian (11). The equations of motion 
follow: For the entropy s, this yields 





(11) 



s = {s, i?}* = -V • Vs, 
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where, as before, we use the 'bar' shorthand for solenoidal quantities, i.e. v = 
V ■ V, and evidently V • v = 0. Note s can be any advected quantity such as the 
concentration of salt. 

Similarly, the dynamical equation for B is obtained 

B = {B, iJ}* = -V • VB + B • Vv = V X (v X B) . 

The equation for v is slightly more complicated, viz. 

V = {v, H}^ = . [(V X v) X v] + P • [(V X B) X B]. (12) 

In particular, the property that V • 7^ = implies that V • v = 0, which is 
consistent with the constraint $2- Wc notice that the first term in Eq. (12) was 
obtained in Rcf. [13]. Since v « v (weak equality with the constraint $2)) the 
equations for v and B becomes 

V = -V • Vv - VPc + (V X B) X B, 
B = V X (v X B), 

where the pressure-like term is given by 
,,2 

Pc := -y - A-^V ■ ((V X v) X v) + A-^V • ((V x B) x B) . 

Given the equation for the pressure, Pc is not necessarily positive. Lastly, we 
point out that there is no equation for the mass density p, since it has been 
eliminated altogether from the theory. 

The equations obtained above correspond to the traditional equations for 
incompressible MHD. It should be noted that V ■ v = is no longer a constraint 
on the flow since it is a conserved quantity. Actually, it is more than a con- 
served quantity since it is a Casimir invariant. If one choses an initial condition 
satisfying V-v 7^ 0, then this quantity will remain constant under the dynamics. 

4. Comparisons between brackets for incompressible fluids 

We focus now on ordinary fluids and in particular we consider different for- 
mulations for describing the motion of an ideal incompressible fluid. In his 
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famous treatise Lagrange [1] provides descriptions of both incompressible and 
compressible ideal fluids. Lagrange uses what is now generally called the La- 
grangian variable description, whereby the dynamics of fluid elements, points, 
are treated in a spatial domain, and he constructs the Lagrangian for this 
inflnite-dimensional system. If we let q denote the position of a fluid element 
labeled by a, that lies in a domain U CM.^ occupied by the fluid, then q: U 
at each time, or q(a, f). To describe incompressible fluids, Lagrange adds the 
constraint det|9q/9a| = 1. It naturally leads to what is now referred to as 
the volume preserving diffeomorphism description of the incompressible fluid. 
A formal description of this was introduced in Refs. [22, 23] for the Euler equa- 
tions of an incompressible fluid. It was based on the fact that diffeomorphisms 
form an infinite parameter Lie group with a Lie algebra given by the commuta- 
tor of vector fields. If V denotes vector fields of K^, then the commutator (Lie 
bracket) 



is again a vector field for any V, W e V, and it is an elementary exercise 
in vector calculus to show the Jacobi identity, [U, [V, W]i,]l+ 0= for all 
U, V, W G V, where denotes the two other terms obtained by cyclic permu- 
tation of (U, V, W). If one restricts T> to contain only divergence- free vector 
fields, V := {V G I>|V • V = 0}, then 'D c I> is a Lie subalgebra, as seen by 
another elementary vector calculation that assures closure: V • [V, Wjj;, = if 
V, W e 'D. From the Lie bracket (13) one can construct the Lie-Poisson bracket 



which indeed satisfies the Jacobi identity for all functionals of v, it being of the 
Lie-Poisson form. However, combined with the Hamiltonian H = it 
does not yield the correct equations of motion for incompressible fluid mechanics 
since V • v is not conserved by the flow. 



[V, W]l = (W • V)V - (V • V)W, 



(13) 



(cf., e.g., Refs. [4, 5]) 




(14) 
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Another bracket for incompressible fluids was proposed in Ref. [24]: 

{F,G}o = jd^xui- [(V X F„) X (V X G^)] . 

With this bracket and the Hamiltonian H = where w = V x v, the 

equation V x H^, = v = leads to 

— = V X (v X u)), 

which is the correct equation of motion for the vorticity in both compressible and 
incompressible barotropic fluids. However, two issues should be noted: (i) this 
bracket does not satisfy the Jacobi identity for functional deflned on arbitrary 
vector flelds oj. This is easily seen by the following counter example: 

Fi = ^ jd^xuj -xy^ , F2^^ jd^xu;-yz^ , , F3 = jd^x u: ■ z x , 

which yields, 

{Fi , {F2, F3}o}o+ 0=-Jd^xuj- S7{yz) + , (15) 

and (ii) it is not stated how the constraint V • v = is to be applied, and 
indeed the procedure of this paper also gives the correct equation of motion for 
compressible barotropic fluids. 

With regards to (i), if one considers vector fields that satisfy w = V x v, 
then Eq. (15) gives zero (see Appendix D). Thus, one might attempt to restrict 
the space of functional on which this bracket is defined in order to get a Lie 
algebra realization on such functional. Yet, since u> = V x v, it seems natural 
just to use V as a variable to enforce the constraint V • a; = 0. This leads to a 
bracket similar to Eq. (14), namely {F, G}, = / (S'xv ■ [Fy, Gv]», where 

[V, W], := V X (V X W) = [V, W]l + V(V • W) - W(V • V). 

This bracket is not of Lie-Poisson type since [•,•], does not satisfy the Jacobi 
identity, as can be seen from the counterexample (Vi, V2, V3) = {xyk,yy,z), 
giving [[Vi, V2],, V3],+ 0= yz. Thus {F, G}, has to be discarded even though 
it has the interesting property {V • v, F}, = for all observables F. 
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This returns us to issue (ii) above about enforcing V ■ v = 0. The failure 
of Eq. (14) to give the correct equations of motion can be traced to the use of 

V X F„ = -Fvj which cannot be true for all functional because = V • V x F^^, = 

V • Fv ^ 0. Note, even if F[v] is defined on divergence-free vector fields, it does 
not follow that V • =0. This suggests introducing 

V X F„ = Fv + T , 

where T is chosen to enforce the constraint, i.e. 

VxF„=Fv-VA-iV-Fv = P-Fv. (16) 

Now, inserting Eq. (16) into Eq. (14), we obtain the Dirac bracket of Sec. 3. 
So the correct Poisson bracket for incompressible fluids can be constructed as a 
Lie-Poisson bracket, from a projection of the Lie bracket [•, -ji, as follows: 

[V,W]p := [P-V,r-W]L, 

where we notice an important property for verifying the Jacobi identity is V ■ 
["P • V, P ■ W]z, = [-P • V, -P ■ W]z, (cf. Appendix C). Previously the need for the 
projection for the incompressible fluid was observed in Refs. [5, 21]. However, 
in light of our work, when projection is handled appropriately, this amounts to 
the Dirac bracket construction of Ref. [13], which we here generalized. 

In closing this section we make a few more remarks. In the two-dimensional 
formulations of Refs. [25, 2, 26] there is no issue with projection: unlike {•, ■}o 
and {•,•}., the bracket given there satisfies the Jacobi identity for all functional 
of the scalar vorticity. Also, in the compressible formulation the density is 
added as a dynamical variable (cf. the first term of Eq. (9) of Ref. [16]) and the 
variations with respect to density in the Jacobi identity compensate the failure 
of Jacobi for the second term alone (see footnotes 10 and 12 of Ref. [16]). Lastly 
we point out that care must be taken when inserting projections on functional 
derivatives into Poisson brackets, for the resulting Poisson bracket may not 
satisfy the Jacobi identity (cf. Appendix B and Appendix C). 
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5. Conclusions 

Here we have generalized the Dirac bracket approach of Ref. [13] by including 
entropy advection. This produces a Hamiltonian description of an important 
missing piece of the dynamics of incompressible fluids, viz. that of density ad- 
vection. Recall, V • V = does not imply constant p, but that p be advected. If 
one chooses V • v as the primary constraint then one does not obtain a bracket 
for an advected density. Thus, it would appear that density advection cannot be 
produced by the Dirac bracket construction. However, having done the calcula- 
tion with entropy advection we observe that p drops out of the picture and we 
obtain a bracket that describes advection of a quantity s by a solenoidal velocity 
field. Thus, if one just reinterprets s as p we obtain the missing dynamics of 
density advection. 

Performing the Dirac construction for MHD with density advection showed 
us that this approach is equivalent to direct projection of the MHD bracket 
of Refs. [16, 2] to solenoidal vector fields. It is now evident how to construct 
brackets for a variety of incompressible models. If it is Lie-Poisson, then one 
can proceed as in Appendix C and if is not, then one can step through the 
Dirac bracket construction. In fact, the Dirac construction is more general and 
can be used to enforce any compatible constraints, as seen, e.g., in Appendix 
A. Evidently, Dirac brackets provide a powerful tool that extends well beyond 
the results of this paper. 
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Appendix A. Generalization to ideal MHD with nonuniform back- 
ground density 

Suppose the density is constant but nonuniform, which might, e.g., be an 
imposed stratification caused by gravity. It is interesting to see where such 
an assumption leads when one follows the Dirac construction. To this end we 
assume $i(x) = p — po(x) = 0, where po is the time-independent background 
density. Proceeding as in Sec. 3, because {4>i(x), <l>i(x')} = we obtain the 
secondary constraint that has the form 

$2(x) = V-(po(x)v) =0. (A.l) 

Although Eq. (A.l) is vahd for compressible equilibria, to justify such a con- 
straint on physical grounds would require a mechanism for maintaining the 
constraint or a time scale argument of some sort. We will not pursue this here. 
From the Poisson bracket (6), we compute the elements Cq/3(x, x') as 

Cn(x,x') = 0, 
Ci2(x,x') = ^(5(x-x'), 

C2i(x,x') = -ASix-x'), 

C22(x,x') = V • {plp-\V X v) X V5(x - x')) , 

where A is the symmetric operator Af = V ■ (poV/). Provided A is invertible, 
we obtain the coefficients C~^(x, x') as 

Cfi^(x, x') = A-^V ■ {plp-\V X v) X Vyt-M(x - x')) , 
C^^\x,^') = -A-'5{x-x'), 
C^,\xy) = A-'S{x-x'), 
(72-/(x,xO = 0. 
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Given the following expressions 
{$i(x),G} = -V-C?v, 

{$2(x), G} = -AGp - V • {pQp-\V X v) X Gv) + V • {p^p-^VsGs) 
-V • (poP"MVGb] • B) + V • (pop-'V • [BGb]) , 

we deduce various contributions to the Dirac bracket (7): 

j j d^xd3a;'{F,$i(x)}GiY(x,x'){$i(x'),G} ^- j d^xV ■ F^A'^V ■ {plp-\V x v) x VA'^V-G^), 

yy^d3xdV{F,$i(x)}Gf2i(x,x'){$2(x'),G} = j (fxV ■ Fv {Gp + A'^V ■ (poP"^(V x v) x Gv) 

-^-^V ■ (pop-'VsG«) + ^-^V • (pop-MVGs] • B) - ^"^V ■ (pop-'V • [BGb])) , 

jj rf3a;dV{J^,$2(x)}G2-/(x,x'){$i(x'),G} = - y d^a; (J^^ + V ■ (pop-'(V x v) x F^) 

-^-^V • {pop-^VsFs)+A-^V ■ (poP"MVFb] • B) - A'^V ■ (poP'^V ■ [BFb])) V • G^ 

The Dirac bracket now reads 

{F, G}* = j d^x (p-i(V X v) • (Fv x Gv) - p-'Vs ■ (f.G^ - F^G,) 

- (p-^Fv • [VGb] - p-'Gv • [VFb]) ■ B - B • ([V (p-^Fv)] • Gb - [V (p-^Gv)] • i^)) 

where =Va-F^ = F^- PoV (^"^ V • Fv) . Observe that V • FV = for these 
equations, as was the case for the incompressible MHD. We also notice that the 
Dirac bracket has the same form as that for incompressible MHD, with the only 
difference being divergence- free functional derivatives F^ replacing F^. 

In the same way as in Sec. 3, one term in the Hamiltonian corresponds to 
a Casimir invariant. More precisely, from the property that V • Fv = for any 
observable F, it is shown that 



C[s] = J d'xpfip, 



is a family of Casimir invariants, where / is any function of s and p. Therefore 
the Hamiltonian is 

H=l [d^x (pov^ + B^) , 



2 , 

and the internal energy U plays no role in the dynamics, just as was the case 
for ideal incompressible MHD. 
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The two dynamical equations for s and B are similar than the ones for 
incompressible MHD, and are given by 

s = —V • Vs and B = V x (v x B) , 

since V = V — pop~^'V{A~^\7 ■ (pv)) « v with the secondary constraint $2- The 
dynamical equation for v becomes 

V = -V • Vv + Pq^{V X B) X B - VWc, 

where the Bernoulli-like term Wc is given by 

Wc = -— - A-'^V ■ {po{V X v) X v) + ^-^V • ((V X B) X B). 

Again we notice that V • (pov) is conserved by the flow since it is a Casimir 
invariant. 

It should be noted that the second constraint $2 above had a constant 
background density. Another choice would be to use the constraint $2 with p 
replacing po, i.e. use the following set of constraints: 

$i(x) = p-po(x) and $2(x) = V • (pv). 

Proceeding as above, the definition of the operator A naturally becomes Af = 
V ■ (pVf), and the expression for the Dirac bracket obtained is identical to 
Eq. (A.2) with := - pV(^-^V • F^), which still satisfles V ■ = 0. 

Appendix B. Hamiltonian-Dirac Vector Fields 

Let Z denote a phase space manifold that is a symplectic or Poisson manifold, 
and is thus equipped with a bracket operation { ■ , • } : C°°{Z) x C°°{Z) — > M. 
We suppose the bracket satisfies the usual Lie enveloping algebra properties and 
can thus be written in coordinates as 

for functions f,g & C°°{Z), i.e. f,g: Z — >• K. Note the bracket above is a 
generic Poisson bracket and may have any form or degeneracy. Only the Lie 
algebra properties are required. 
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We impose an even number of constraints G C°°{Z), a = l,. . ., 2m, and 
wish to project Hamiltonian vector fields on Z, elements of X{Z), to Hamilto- 
nian vector fields that are tangent to a submanifold Ai := fla^a, X{M). 

As is well-known elements of X{Z) are linear operators, in particular, the 
element generated by / G C°° (Z) has the form 

r - If \- rabdf d 
Lf--U,-}-J 

and the commutator of two such elements satisfies [Lf,Lg] = —L^f gy. Thus 
there is an isomorphism between the Lie algebra of such linear operators and 
Poisson brackets. We wish to maintain this structure for Hamiltonian vector 
fields projected onto X{A4). 

To project a Cartesian vector into a surface defined by (f) =constant, one 
uses the normal Wcj) to construct the following projection operator: 

where I is the identity. Evidently P • = 0. Essentially this same idea 
occurs in infinite dimensions in the context of Hilbert spaces and is efficacious 
for application in quantum mechanics. However, the problem at hand differs 
from these cases in that we are interested in Hamiltonian vector fields (finite or 
infinite) and our manifold is symplectic with no intrinsic notion of metric. Thus, 
if we are to proceed without adding additional structure, we must construct a 
projection operator using only the functions and cosymplectic form, J. With 
Eq. (B.l) as a guide we write 

where K^p is chosen so that Hamiltonian vector fields generated by any of the 
$a are projected out, i.e. V ■ Z/$^ = for all a. Now it is desired to find such a 
in terms of the and J alone. Fortunately, a direct calculation reveals 
that the desired quantity is given by = {^a,^p}~^ ■ Thus we have achieved 
our goal if this inverse exists. Assuming this is the case we obtain the following 



16 



Hamiltonian projection operator: 



n^n- 



Evidently 



(B.2) 

and = for all a. Also, an elementary calculation reveals the = V, as 
expected for a projection operator. 

It remains to show that the set of projected vector fields of the form Aj = 
V ■ Lf are Hamiltonian on the constraint submanifold: 

[A/,A,] = -A{/,,},, (B.3) 

for some well-defined Poisson bracket {/, 5}*. As the notation suggests this 
turns out to bo the Dirac bracket. 

It is evident from Eq. (B.2) that h"-^ — — {/, •}*. Because Eq. (B.3) is 
satisfied for a generic Poisson bracket, then it must be true for the Dirac bracket 
as well. To see that it is true for a generic Poisson bracket we write 

= -ir'dbr^djd.g + j'^r'dhidjdsg) 

where db := d/dz^ and db operates only on the term immediately to its right un- 
less parenthesis are included. In obtaining the second equality, second derivative 
terms canceled in the usual way, and in obtaining the third equality, antisym- 
metry, the Jacobi identity, and relabeling were used. 

All of the above can be formally extended to infinite dimensions (see, e.g., 
Ref. [2]) by replacing partial derivatives by functional derivatives, sums by in- 
tegrals, and matrix multiplication by operator action. 
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Appendix C. Projections and Poisson brackets 

Consider the general Poisson bracket, 

where J is a cosymplectic operator (generally dependent on x(/x)) that ensures 
this bracket satisfies the Jacobi identity. Now suppose V is some projection 
operator, and consider 

The bracket (C.l) does not in general satisfy the Jacobi identity. However, if V 
is independent of x it may. 

If the bracket is Lie-Poisson, then projection onto subalgebras always pro- 
duces brackets that satisfy the Jacobi identity. Consider the Lie-Poisson bracket 



{F,G} = Jd^,^^]S^-^ = {x,[F^,G^])^ 



In this construction & g where g is a Lie algebra and hence is a vector. 
Suppose V : Q ^ t, where fi is a vector subspace of q. Then, (i) the bracket 

{F,G}p = {x,['PF^,VG^]), (C.2) 

is defined, and (ii) it satisfies the Jacobi identity for arbitrary functional of Xj 
provided VlVF^jVG^] = [PF^.VG-^, which is the case if 6 is a subalgebra of 
0. This follows from the general Jacobi identity theorem proven in Ref. [2] or 
more immediately from the fact that Eq. (C.2) is a Lie-Poisson bracket for 6. 

Appendix D. Direct proof of Jacobi identity 

We consider the bracket 

{F,G] = J d^xu:- [(V X F^) x (V x G^)] (D.l) 

In order to prove the Jacobi identity for this kind of bracket, one only needs to 
consider the explicit dependence of the bracket on u> when taking the functional 
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derivative 6{F,G}/du}. In what follows, we let f := V x The functional 
derivative of {F,G} with respect to u> contains three terms: one that comes 
from the explicit dependence of the bracket on w and two other terms that 
are the second order functional derivatives of F and G. It has been shown in 
Ref. [2] that the only important term comes from the explicit dependence on 
the variables, i.e. on a>. So from {F, G}^, = f x g, we get 

{F,{G,H}} = J d^xu: ■ (f X V X (g X h)). 

Since V • f = 0, this becomes 

{F, {G, if }} = y w ■ [f X (h ■ V)g - f X (g • V)h] . (D.2) 

If V ■ u> = 0, there exists a vector v such that u> = V x v. By symmetry of the 
operator Vx, we obtain terms like V x [f x (h • V)g — f x (g • V)h], which are 
transformed by the identity 

V X [f x (h • V)g] = fV • [(h • V)g] + [(h • V)g • V]f - (f • V)[(h • V)g]. 

Since V • f = (for all observable F), the divergence terms in Eq. (D.2) vanish, 
i.e. V • [(h • V)g] = V • [(g • V)h]. In addition, from the following identity 
(f • V)[(h • V)g] = [(f ■ V)h • V]g + fihjdidjg, we have 

Vx[fx(h.V)g-fx(g.V)h] = (h,g,f)-(f,h,g)-(g,h,f)+(f,g,h)-/i/i,aia,g+/iff,-5i5,h, 

where (f , g, h) := [(f • V)g • V]h. By adding the cycUc permutations of F, G 
and H, we obtain 

V X [f X (h • V)g - f x (g • V)h]+ 0= 0. 

As a consequence, the bracket (D.l) satisfies the Jacobi identity if V • w = 0. 

From the bracket (D.l) with V • w = 0, we perform the following change of 
variables: a> = V x v. This change of variables depends on a gauge since v + W(f) 
gives the same value for us. For instance, if we choose V • v = 0, we obtain 

V xF^^F^- VA-iy • F^. 

Thus, we end up with the Poisson bracket obtained with Dirac's procedure for 
constrained Hamiltonian systems. 



19 



References 

[1] J.L. Lagrange, Mecanique Analytique, English Title: Analytical Mechan- 
ics, translated and edited by A. Boissonnade and V.N. Vagliente (Kluwer 
Academic, Imprint Dordrecht, Boston, Mass. 1997). 

[2] P.J. Morrison, in Mathematical Methods in Hydrodynamics and Integrability 
in Related Dynamical Systems, AIP Conference Proceedings 88, edited by 
M. Tabor and Y. Treve (AIP, New York, 1982), p. 13. 

[3] R. Salmon, Annu. Rev. Fluid Mech. 20, 225 (1988). 

[4] P.J. Morrison, Rev. Mod. Phys. 70, 467 (1998). 

[5] J.E. Marsden and T.R. Ratiu, Introduction to Mechanics and Symmetry 
(Springer- Verlag, Berlin, 2002). 

[6] P.J. Morrison, Phys. Plasmas 12, 058102 (2005). 

[7] P.J. Morrison, "Hamiltonian Fluid Mechanics", in Encyclopedia of Mathe- 
matical Physics, vol. 2, (Elsevier, Amsterdam, 2006) p. 593. 

[8] R. Salmon, J. Fluid Mech. 196, 345 (1988). 

[9] J. Vanneste and O. Bokhove, Physica D 164, 152 (2002). 

[10] P.J. Morrison, N.R. Lebovitz and J.A. Biello, Ann. Phys. 324, 1747 (2009). 

[11] C. Chandre, E. Tassi and P.J. Morrison, Phys. Plasmas 17, 042307 (2010). 

[12] G. Flierl and P.J. Morrison, Physica D 240, 212 (2011). 

[13] S. Nguyen and L.A. Turski, Physica A 272, 48 (1999). 

[14] S. Nguyen and L.A. Turski, Physica A 290, 431 (2001). 

[15] S. Nguyen and L.A. Turski, Physica A 388, 91 (2009). 

[16] P.J. Morrison and J.M. Greene, Phys. Rev. Lett. 45 , 790 (1980); ibid. 48, 
569 (1982). 

20 



[17] P.A.M. Dirac, Can. J. Math. 2, 129 (1950). 

[18] E.C.G. Sudarshan and N. Mukunda, Classical Dynamics: A Modem Per- 
spective (John Wiley & Sons, New York, 1974). 

[19] A. Hanson, T. Rcggc and C. Tcitclboim, Constrained Hamiltonian Systems 
(Accademia Nazionale dei Lincei, Roma, 1976). 

[20] K. Sundermeyer, Constrained Dynamics (Springer- Verlag, Berlin, 1982). 

[21] V.E. Zakharov and E.A. Kuznetsov, Physics-Uspekhi 40, 1087 (1997). 

[22] V.I. Arnold, Ann. Inst. Four. 16, 319 (1966). 

[23] V.I. Arnold, Usp. Mat. Nauk. [Sov. Math. Usp] 24, 225 (1969). 

[24] E.A. Kuznetsov and A.V. Mikhailov, Phys. Lett. 77A, 37 (1980). 

[25] P.J. Morrison, "Hamiltonian field description of two-dimensional vortex 
fluids an guiding center plasmas", Princeton University Plasma Physics 
Laboratory Report, PPPL-1783 (1981). Available as American Institute of 
Physics Document No. PAPS-PFBPE-04-771-24, AIP AuxiUary Publica- 
tion Service, 335 East 45th Street, New York, NY 10017. 

[26] P.J. Olver, J. Math. Anal. Appl. 89, 233 (1982). 



21 



